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Abstract 

We show that the minimal R-parity breaking model characterized by an effective 
bilinear violation of R-parity in the superpotential is consistent with minimal N=1 
supergravity unification with radiative breaking of the electroweak symmetry and 
universal scalar and gaugino masses. This one-parameter extension of the MSSM- 
SUGRA model provides therefore the simplest reference model for the breaking of 
R-parity and constitutes a consistent truncation of the complete dynamical mod¬ 
els with spontaneous R-parity breaking proposed previously. We comment on the 
lowest-lying CP-even Higgs boson mass and discuss its minimal N=1 supergravity 
limit, determine the ranges of tan/3 and bottom quark Yukawa couplings allowed in 
the model, as well as the relation between the tau neutrino mass and the bilinear 
R-parity violating parameter. 





1 Introduction 


Supersymmetry apart from being attractive from the point of view of providing a solution 
to the hierarchy problem and the unihcation of the gauge couplings, provides an elegant 
mechanism for the breaking of the electroweak symmetry via radiative corrections [Q. So 
far most attention to the study of supersymmetric phenomenology has been made in the 
framework of the Minimal Supersymmetric Standard Model (MSSM) [Q with conserved 
R-parity R-parity is a discrete symmetry assigned as Rp = (—l)( 3 -B+i+ 25 )^ where L is 
the lepton number, B is the baryon number and S is the spin of the state. If R-parity is 
conserved all supersymmetric particles must always be pair-produced, while the lightest 
of them must be stable. Whether or not supersymmetry is realized with a conserved 
R-parity is an open dynamical question, sensitive to physics at a more fundamental scale. 

The study of alternative supersymmetric scenarios where the effective low energy 
theory violates R-parity has received recently a lot of attention i,i0- As is well- 
known, the simplest supersymmetric extension of the Standard Model violates R-parity 
through a set of cubic superpotential terms involving a very large number of arbitrary 
Yukawa couplings. Although highly constrained by proton stability, one cannot exclude 
that a large number of such scenarios could be viable. Nevertheless their systematic study 
at a phenomenological level is hardly possible, due to the large number of parameters 
(almost fifty) characterizing these models, in addition to those of the MSSM. 

As other fundamental symmetries, it could well be that R-parity is a symmetry 
at the Lagrangean level but is broken by the ground state. Such scenarios provide a 
very systematic way to include R parity violating effects, automatically consistent with 
low energy baryon number conservation. They have many added virtues, such as the 
possibility of having a dynamical origin for the breaking of R-parity, through radiative 
corrections, similar to the electroweak symmetry ^j. 

In this paper we focus on the simplest truncated version of such a model, in which the 
violation of R-parity is effectively parametrized by a bilinear superpotential term CiL^H^. 
In this effective truncated model the superfield content is exactly the standard one of the 
MSSM. In this case there is no physical Goldstone boson, the Majoron, associated to the 
spontaneous breaking of R-parity. Formulated at the weak scale, the model contains only 
two new parameters in addition to those of the MSSM. Alternatively, the unified version 
of the model, contains exactly a single additional parameter when compared to the unihed 
version of the MSSM, which we will from now on call MSSM-SUGRA. Therefore our model 
is the simplest way to break R-parity and can thus be regarded as a reference model for 
R-parity breaking. In contrast to models with trilinear R-parity breaking couplings, it 
leads to a very restrictive and systematic pattern of R-parity violating interactions. 

Here we show that this simplest truncated version of the R-parity breaking model 
of ref. PI, characterized by a bilinear violation of R-parity in the superpotential, is 
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consistent with minimal N=1 snpergravity models with radiative electroweak symmetry 
breaking and nniversal scalar and gangino masses at the nnihcation scale. In particular, 
we perform a thorough study of the minimization of the scalar boson potential using the 
tadpole method needed for an accurate determination of the Higgs boson mass spectrum. 
We comment on the lowest-lying CP-even Higgs boson mass and discuss its minimal 
N=1 snpergravity limit, determining also the ranges of tan/? and bottom quark Yukawa 
couplings allowed at nnihcation, as well as the relation between the tan neutrino mass 
and the effective bilinear R-parity violating parameter. Our results encourage further 
theoretical work on this and on more complete versions of the model, like that of ref. P] , 
as well as phenomenological studies of the related signals. 


2 The Model 


The supersymmetric Lagrangian is specihed by the superpotential W given by Q 

W = \h^Q'iU,H^2 + hWiDjHl + h%L%H^ - ( 1 ) 


where i,j = 1,2,3 are generation indices, a,b = 1,2 are SU{2) indices, and £ is a com¬ 
pletely antisymmetric 2x2 matrix, with ei2 = 1. The symbol “hat” over each letter 
indicates a superheld, with Qi, Li, Hi, and H 2 being SU{2) doublets with hyper-charges 
—1, —1, and 1 respectively, and U, D, and R being SU{2) singlets with hyper-charges 
— and 2 respectively. The couplings hjj, ho and He are 3x3 Yukawa matrices, and 
and Cj are parameters with units of mass. The hrst four terms in the superpotential are 
common to the MSSM, and the last one is the only i?-parity violating term. From now 
on, we work only with the third generation of quarks and leptons. 


Experimental evidence indicate that supersymmetry must be broken. The actual 
snpergravity mechanism is unknown, but can be parametrized with a set of soft supersym¬ 
metry breaking terms which do not introduce quadratic divergences to the unrenormalized 


theory [10 


Vsoft = m^qQ'^,*q^, + m^u;us + mId*Ds + mIl'^,*l'^, + mIrirs 




-|- XM 2 X 2 X 2 hMiXiXi h.c. (2) 




+ A,KL\R^HI - + B2e^LlHl 


where we are already using a one-generation notation. 

Note that in the effective low-energy snpergravity model the bilinear R-parity vio¬ 
lating term cannot be eliminated by superheld redehnition even though it appears to be 
so at high scales, before electroweak and supersymmetry breaking take place The rea¬ 
son is that the bottom Yukawa coupling, usually neglected in the renormalization group 

* We are using here the notation of refs. ||] and |^. 
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evolution, plays a crucial role in splitting the soft-breaking parameters B and B2 as well 
as the scalar masses and M|, assumed to be equal at the unihcation scale. This can 
be seen explicitly from eq. (|^) and eq. (^7]) as well as eq. (|i9[) and eq. ( ^21) in Appendix 
A. This ensures that R-parity violating effects can not be rotated away by going to a new 
basis □ 0,0, even if the starting RGE boundary conditions for the soft-breaking terms 
are universal. 


It goes without saying that, in a supergravity model where soft-breaking terms are 
not universal at the GUT scale, such as string models, the bilinear violation of R-parity 
is also not removable. However, in this case its effects are not calculable, in contrast to 
our case. The same is true for the case of the most general low-energy supersymmetric 
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The electroweak symmetry is broken when the two Higgs doublets Hi and H 2 , and 
the tau-sneutrino acquire vacuum expectation values (VEVS): 


Hi = 




Ho = 


Ht 


,^[X2 + V2+W2], 




1 + ^,3 + 


T~ 


( 3 ) 


Note that the gauge bosons W and Z acquire masses given by and m| = 

\{g'^ + g'‘^)v'^ 1 where v‘^ = vl + = (246 GeV)^. We introduce the following notation 

in spherical coordinates: 

vi = V sin 6 cos {3 
V 2 = V sin 9 sin (3 

V 3 = vcos9 (4) 


which preserves the MSSM dehnition tan/? = V 2 IV 1 . The angle 6 equal to 7r/2 in the 
MSSM limit. 


The full scalar potential may be written as 




dW 

dzi 


2 

+ Vd + Vsoft + Vrc 


( 5 ) 


where Zi denotes any one of the scalar helds in the theory, Vd are the usual H-terms, 
Vsoft the SUSY soft breaking terms given in eq. (^, and Vrc are the one-loop radiative 
corrections. It is popular to treat radiative corrections with the effective potential. In this 
case, Vrc corresponds to the one-loop contributions to the effective potential. Here we 
prefer to use the diagrammatic method and hnd the minimization conditions by correcting 
to one-loop the tadpole equations. At the level of Ending the minima, the two methods 

t Obviously physics does not depend on the choice of basis Q. In this paper we choose to work with 
the unrotated fields. 
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are equivalent m- Nevertheless, the diagrammatic (tadpole) method has advantages 
with respect to the effective potential when we calculate the one-loop corrected scalar 
masses 


The scalar potential contains linear terms 


Vlinear = tiXl + ^ 2^2 + ^ 3 ^- 


0^,0 


M~R 


( 6 ) 


where 

+ |(^^ + g^)vi{vl - + vl), 

tl = + 63)^2 - B/j^vi + ^ 263^3 - |(^^ + 9’^)v2{vI -vi + vl ), (7) 

^3 = i^L3 + 4)^3 - he3^'l + B2e3V2 + 1(5'^ + 9'‘^)v3{vI -vi + Vs) . 

These i = 1, 2, 3 are the tree level tadpoles, and are equal to zero at the minimum of 

the potential. 


3 Squark Sector and Radiative Corrections 


To hnd the correct electroweak symmetry breaking radiatively, we need to relate param¬ 
eters at the GUT scale with parameters at the weak scale. This means we are promoting 
the parameters in the tree level tadpoles in eq. (0) to running parameters. Therefore, 
in order to hnd the correct minima of the scalar potential it is essential to include the 
one-loop contributions to the tadpoles, otherwise our tadpoles would be extremely scale 
dependent, he., unphysical. 


The main one-loop contributions to the tadpoles come from loops involving top and 
bottom quarks and squarks. Therefore, we need to study the scalar quark sector, and in 
particular, the spectrum and couplings to CP-even neutral scalars. 

The term esL^H^ in the superpotential induce F-terms in the scalar potential, lead¬ 
ing to squark mass terms of the form proportional to 63 . In addition, the non-zero 
value of the vacuum expectation value of the tau-sneutrino generates, from the D-terms, 
squark mass terms of the form ui* and bib*, i = L,R. The new squark mass matrices are: 


Ml 


Mq + ml + + vl) mt{At - 9.V1/V2 + esVs/v2) 

mt{At - HV 1 /V 2 + 63 ^ 3 ^) MIj -h ml + \9'‘^{vl - + vfj 


( 8 ) 


for the top squarks, and 


Ml = 




+ ml 


-U9^ + 




){vl-vl + vl) 


mb{Ab - 9 ,V 2 /vi) 


mb{Ab - 9LV2/V1) 


Ml) -h ml 


y\vi-vi+vi) 


(9) 


for the bottom squarks. The reader can recover the MSSM squark mass matrices by 
taking 63 = ^3 = 0 in the above two equations. The quark masses are related to the quark 
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Yukawa couplings in the same way as in the MSSM: rrit = htV2/V2 and rrib = hf,vi/y/2. 
Nevertheless, the numerical value of the quark Yukawas is higher in comparison with the 
MSSM to compensate with smaller vacuum expectation values 

u - _ grrih 

^ 1 hf) — , ( 10 ) 

y Zmy/SfjSo y 2m\YCpSg 

and this is represented by the term sin 6^ = sg in the denominators in the above equations. 


Squark mass matrices Mf and Mf are diagonalized by two rotation matrices such 


that: 


RiMWi = 


mil 

0 


0 


m 


^2 


RMRl = 


m 


hi 


0 


0 




( 11 ) 


where < mq^ by convention. These rotation matrices play an important role in the 
determination of the scalar couplings to a pair of squarks. 

We introduce the notation for the CP-even neutral scalars = XiiX 2 i^t 1°^ ^ = 
1, 2, 3 respectively. In this way, the Feynman rules of the type S^qq are 


x\bh 








( 12 ) 


as in the MSSM, but with the quark Yukawa couplings given by eq. m- Feynman rules 
of the type S^qq not listed in eq. 0 ) are zero. 

In a similar way, we hnd Feynman rules of the type Sfqq*, i.e., CP-even neutral 
scalars couplings to a pair of squarks. We start with Xi couplings to top squarks: 


Xiii* 

^ Xitt ~ 


’ ^Xiit ~ Ri^ XitiRi ’ 




75 


-^htH 


P 2 




(13) 


and to bottom squarks: 

xibb' 


^^Xihh ’ 


— RiM , 


^y^h= 


Xibb — xibb-^b ’ 

-hl^i + lig'^ + y‘^)vi -^hbAb 


4 Vi< ' 3- 
-^hbAb 


-Hvi + ^g'^vi 


(14) 


These couplings have the same form in the MSSM but, as it was said before, the Yukawa 
couplings are different and given by eq. (^). In addition, vacuum expectation values 
Vi and V 2 are different with respect to the MSSM and given by Vi = 2 mwCf3SQ/g and 
V 2 = 2mwSj3Se/g and again, the deviation from the MSSM is parametrized by the angle 
6 . 

Now we turn to the neutral CP-even Higgs X 2 comes from the second Higgs 
doublet. Its couplings to top squarks are: 


X 2 ^t* 


M xitt 


^xlti ~ Ri^ xiiiRi ’ 

-h^tV 2 + lig^ - y^)v 2 




-^htAt -htV 2 + y^V 2 


(15) 
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and to bottom squarks: 


x>’ 


iM 


^Y°bb — ^b^ X2bb^b ’ 


X°2bb 


xfbb^ ^''^xfbb 
■^hb^i 
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( 16 ) 


Finally, we tnrn to the real part of the tan-sneutrino held, which mixes with Xi and 


X 2 - Its conplings to top squarks are: 


^ iM 


u^tt 
^'mi = 


v^it 


^vMi — vMt^i 1 


-jig^ - 2g'^)v3 --jnhtes 






and to bottom squarks: 


(17) 




^ uPbb ~ 


iM 


u^bb 1 

\{g^ + \g>, ’ 0 
0 y^v. 


^i>Rbb — Rf,^ uRbbRj, ’ 


These couplings p^qq* vanish in the MSSM limit = 0, as it should. 


(18) 


We are now ready to include the effect of the one-loop tadpoles in eq. 0 . The 
hrst step towards the calculation of radiative corrections is the introduction of counter¬ 
terms. All parameters in the Lagrangian are shifted from bare parameters to renormalized 
parameters minus a counter-term: 


A —^ A - 5A 

—>■ 

Vi —^ Vi - 6vi 

A —^ A-SA 
B —^ B-5B 


A g 1 g ) ht, hb, h-T, 

* = 1,2,3, 

A = At, Ah, At, 

B = B, i? 2 . 


(19) 

( 20 ) 


for couplings, masses, vacuum expectation values, trilinear soft parameters, and bilinear 
soft parameters respectively. If we make this shift in the tadpole equations given in 
eq. (0), the tadpole themselves get a counter-term 6ti for i = 1,2,3. Therefore, the 
one-loop tadpole equations are 

ti = tySti + Ti{Q), ^ = 1,2,3, (21) 

where ti are the one-loop renormalized tadpoles and Ti{Q) are the one-loop contribu¬ 
tions to the tadpoles, with Q being the arbitrary mass scale introduced by Dimensional 
Reduction. 


The renormalization scheme we choose to work with is the MS scheme, where by 
definition the tadpole counter-terms are taken such that they cancel the divergent pieces 
of Ti{Q) proportional to A: 

A =-^-hln47r-7E, (22) 

4 — n 
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where A is the regulator of dimensional regularization, n is the number of space-time 
dimensions, and 7 ^; is the Euler’s constant. The MS'-counter-terms chosen in this way 
make the tadpoles hnite. We introduce the notation 

rp(g) =-<5tp + T,(Q), (23) 


for the hnite one-loop contribution to the tadpoles. These hnite one-loop tadpoles depend 
explicitly on the arbitrary scale Q. 

The one-loop tadpoles ti must be scale independent (at least in the one-loop approx¬ 
imation), therefore, the renormalized parameters are promoted to running parameters, 
he., they evolve with the scale Q according to their Renormalization Group Equations 
(RGE). The explicit Q dependence on T^^{Q) is cancelled at one-loop by the implicit 
Q dependence on the parameters of the tree level tadpoles. Renormalized tadpoles must 
be zero at the minimum of the potential ti = 0 , thus the generalization of the tadpole 
equations is 


{m]j^ + - BfiV2 - /iesus + - 

+ g? + el)v2 - Bgvi + R2e3^'3 - |(fi'^ + g''^)v2(vf 
(mig -f el)vs - /icsui -f + |(5'^ + g'‘^)v3{vf - 


vl + (Q) + TF^iQ) = 0 , 

-vl + vD] {Q)+fp{Q) = t), 

Vl + (Q) + = 0 . (24) 


and these are the minimization condition we impose 0. We choose to work at the scale 
Q = ruz- The RGE’s for each parameter are given in the Appendix A, and the boundary 
condition at the GUT scale are described later. 


Now we hnd the one-loop contributions to the tadpoles. Quarks contribute to Xi 
and X 2 one-loop tadpoles only. On the contrary, squarks contribute to all three tadpoles. 
Using the notation for the Feynman rules introduced in the previous section, the quark 
and squark one-loop contribution to the tadpoles can be written as: 


^x? 


'^xl 


Tr, 


tbtb 


tbtb 


tbtb 


Nr. 


167r2 




Nr 


167r2 




+ 


+ 


Ncgml 

Sir'^mwCisSe 

Ncgrrij 

Sn'^mwSisSe 


Mml) 

Ao(m2) 


Nr 




167r2 

where Aq is the hrst Veltman’s function dehned by 


(25) 


m 


Ao{m ) = m (A - In — -h 1 ) 


(26) 


The hnite tadpoles T^^{Q) are found simply by setting A = 0 in the previous expressions. 

^ To see the effect one-loop tadpoles have on the determination of MSSM-SUGRA parameters, see 
ref. IT^ 
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4 Unification 


We now discuss the corresponding boundary conditions at unification. We assume that 
at the unification scale the model is characterized by one universal soft supersymmetry¬ 
breaking mass mo for all the scalars (the gravitino mass), and a universal gaugino mass 
Mi/ 2 - Moreover we assume that there is a single trilinear soft breaking scalar mass 
parameter A and that the bilinear soft breaking parameter B is related to A through 
B = A — 1. In other words we make the standard minimal supergravity assumptions: 


cr 

III 

1^ 

(27) 

B = B2 = A -1 , 

(28) 

= Ml = Ml = ml , 

(29) 

Ml = Ml = Ml = ml, 

(30) 

= A /2 “ — Ml ^2 

(31) 


at Q = Mcut- At energies below Mgut these conditions do not hold, due to the renor¬ 
malization group evolution from the unification scale down to the relevant scale. 


In order to determine the values of the Yukawa couplings and of the soft breaking 
scalar masses at low energies we first run the RGB’s from the unification scale Mgut ~ 
10^® GeV down to the weak scale. In doing this we randomly give values at the unification 
scale for the parameters of the theory. The range of variation of the MSSM-SUGRA 
parameters at the unification scale is as follows 


10-2 

< 


< 

1 


10“^ 

< 


< 

1 


-3 

< 

A/mo 

< 

3 

(32) 

0 

< 

Agut/^o 

< 

10 


0 

< 

Mi/ 2/mo 

< 

5 


' es is 






10- 

< 

2 / 2 

< 10 


(33) 


and the value of is defined in such a way that we get the r mass correctly. 

After running the RGB we have a complete set of parameters, Yukawa couplings and 
soft-breaking masses mj{RGE) to study the minimization. 


Similar to what happens in the MSSM-SUGRA (see Appendix B) the number of 
independent parameters of this model is actually less than given above, as one must take 
into account the W mass constraint as well as the minimization conditions. In the end 
there is a single new parameter characterizing our model, namely € 3 . 



5 Results and Phenomenology 


The main parameters characterizing electroweak breaking are the SU(2) doublet VEVs 
Ui, V 2 and V 3 . In our model these are obtained as explained in the Appendix B. Basically 
we assign random values for the top and bottom quark Yukawa couplings ht and h}, at the 
GUT scale and evolve them down to the weak scale through the Renormalization Group 
Equations, given in Appendix A. Using the measured top and bottom quark masses we 
determine the corresponding running masses at the weak-scale. Gombining this with the 
values of ht and hb at the weak-scale, obtained through the use of the RGE’s, we calculate 
the standard MSSM VEVS Vi and V 2 - The third VEV U 3 , which breaks R-parity, is 
determined through the W mass formula. The resulting VEVs may not be consistent 
with the minimization conditions. In Appendix B we present a procedure to ensure a 
consistent solution. Note that due to the contribution of U 3 to the intermediate gauge 
boson masses, v\ + v\ is smaller than in the MSSM. The first check of we can do to 
verify the consistency of the model is to study the allowed values of the lightest GP-even 
Higgs boson mass mb, as a function of the third VEV U 3 . This is displayed in Fig. (|^) 
The unrotated neutral GP-even Higgs bosons Xi cind X 2 with the real part of the tau 
sneutrino These are the GP-even scalars S'?, i = 1, 2, 3, introduced in section 3. The 
mass matrix can be written as 


M|o — ^'s°,MSSM + ^ 


50,63 


+ M' 


S 0 ,RC 


(34) 


where M|o mssm MSSM mass matrix given by 


S°,MSSM 


(35) 


+ zhzn 0 

-Bfi - \glv1V2 0 

0 0 rnl + yUvf-vl)\ 

where we have defined = g^+g'"^. As expected, this mass matrix has no mixing between 
the Higgs and stau sectors. The extra terms that appear in our cs-model are 


M go — 


63 


/ies 


U3 


+ Igl^iv^ 


(36) 


0 

0 -^263f ^263 - \glv2V2, 

+ lglviV3 3263 - \glv2V3 

which introduce a Higgs-Stau mixing. Finally, in M|o introduce the largest term 

in the one-loop radiative corrections, he., the term proportional to 


A- — 7 ?^ 

^3 ^ 8 ^^^3 


^S 0 ,RC 


0 0 0 
0 At 0 
0 0 0 


At = 


3g‘^mj 




In 


2 2 
m- m~ 

tl 12 . 
4 ’ 

mf 


(37) 


This formula gives results good in first approximation, nevertheless, already in the MSSM 
can give wrong results in certain regions of parameter space [O , and should be improved. 


9 













/ (MSSM) 






^140 

0 

o 

E 

120 

100 


80 


60 

1 10 10 ^ 

tan p 

Figure 2: Lightest CP-even Higgs boson mass nih versus tan/? 

As one can see in Fig. (|l|), in the limit Us —> 0 our model reproduces exactly the 
expected minimal SUGRA limit for the lightest CP-even Higgs boson mass. Another 
view of the Higgs boson mass spectrum allowed in our model is obtained by plotting 
rrih as a function of tan/?, as illustrated in Fig. (Q). One sees that all values of tan/? 
in the range 2 to 60 or so are possible in our model. As in the MSSM-SUGRA, tan/? 
smaller than 2 are not possible because the top Yukawa coupling diverges as we approach 
the unihcation scale. This is related to the fact that in that region we are close to the 
infrared quasi-hxed point. Note that the range of tan/? values obtained in our model is 
consistent with the unihcation hypothesis for a large range of the bottom quark Yukawa 
coupling at unihcation, as illustrated in Fig. (^. 

Another important feature of our broken R-parity model is that the tan neutrino 
v-r acquires a mass, due to the fact that es and are nonzero. Consider the basis 
= (—zAi, —zA 2 , i^r), where Ai is the U{1) gaugino introduced in eq. (||), A 2 

is the neutral SU{2) gaugino, iL*, z = 1,2 are the neutral Higgsinos, and Ur is the SM 
tan neutrino. In this base, the mass terms in the Lagrangian for the neutralino-neutrino 
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Figure 3: tan/? versus bottom quark Yukawa coupling at unification 
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sector are 



r 

= 


-b h.c. 


where the mass matrix is 

0 

■ Ml 

0 


\9'v2 

-yvz' 


0 

M 2 

\9Vi 

-\9V2 

\9v^ 

Mjv = 

-\g'vi 

\gvi 

0 

-h 

0 


\g'v2 

-\gv2 

-h 

0 

es 


--\g'v3 

\9V:i 

0 

es 

0 . 


( 38 ) 


(39) 


The only new terms appear in the mixing between neutralinos and tau-neutrino. This 
mixing is proportional to 63 and ^ 3 . They lead to a non-zero Majorana v-r mass, which 
depends quadratically on the lepton-number-violating parameters 63 and ^ 3 . Thus R- 
parity violation in this model is the origin of neutrino mass. In Fig. (^) we display 
the allowed values of rriy^ (in the tree level approximation) as a function of an effective 
parameter ^ dehned as .^ = ( 63^1 -f Notice that rriy^ values can cover a very wide 

range, from eV to values in the MeV range, comparable to the present LEP limit . 
The latter places a limit on the value of € 3 . Note that the values of ^3 and 63 can be 
rather large [see, for example. Fig. (Q)]. 


6 Discussion and Conclusions 


Here we have shown that this simplest truncated version of the R-parity breaking model of 
ref. , characterized by a bilinear violation of R-parity in the superpotential, is consistent 
with minimal N=1 supergravity models with radiative electroweak symmetry breaking 
and universal scalar and gaugino masses at the unihcation scale. We have performed 
a thorough study of the minimization of the scalar boson potential of the model, using 
the tadpole method. We have determined the lowest-lying CP-even Higgs boson mass 
spectrum. We have discussed how the minimal N=1 supergravity limit of this theory is 
obtained and verihed that it works as expected. We have determined also the ranges of 
tan (3 and bottom quark Yukawa couplings allowed at unihcation, as well as the relation 
between the tau neutrino mass and the effective bilinear R-parity violating parameter. 
Our results should encourage further theoretical work on this model, as well as more 
complete versions of the model, like that of ref. |p. Phenomenological studies of the 
related signals should also be desirable, given the fact that the production and decay 
patterns of Higgs bosons and supersymmetric particles in this model are substantially 
different that expected in the MSSM or MSSM-SUGRA. For example, Higgs bosons may 
have sizeable R-parity violating decays |^. Similarly, sneutrinos and staus can be the 


§More complete forms of this matrix have been given in many places. See, e.g. ref. [ p^ 

^This combination appears in treating the neutral fermion mass matrix in the seesaw approximation. 
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LSP and can have nnsnppressed decays into standard model states, thns violating R- 
parity. Finally, chargino and nentralino prodnction can lead to totally different signals 
as, for example, the lightest nentralino can decay 0. These featnres conld play an 
important role in designing strategies for searching for snpersymmetric particles at fntnre 
accelerators. For example, R-Parity will give rise to enhanced lepton mnltiplicities in 
Glnino Cascade Decays at LHC . 
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Appendix A: The Renormalization Group Equations 


Here we will give the renormalization group equations for the model described by the 
superpotential in eq. (|^), but including only the third generation, and by the soft super- 
symmetry breaking terms given in eq. (§). First we write the equations for the yukawa 
couplings of the trilinear terms: 

= hu {qHIj + - 

= ho + hi 

= K {K + - 3gl - 3gl) (42) 

now the corresponding cubic soft supersymmetry breaking parameters 

= QhljAu + HIjAd + —g^M^ -|- 3glM2 -|- —glMi (43) 

= QhljAo -|- hljAjj -|- HIAt- + —glM^ + 3glM2 + -glMi (44) 

HA 

Svr^-^ = 4:hlAr + 3hljAD 3glM2 + 3glMi (45) 

For the soft supersymmetry breaking mass parameters we have 


16 n ^2 13 2\ 

~^93 ~ ^92 

( 40 ) 

16 2 Q 2 2\ 

- y^3 -3^2 - g^ij 

( 41 ) 


Svr^ 


,dM^ 

dt 


— + ^u) + + Mq + Mh Ajj) 


~glMl - 3glMl - ^glMf + ^ glS 


+ A^) - - ^-^glMl - \ glS 

Stt^ —7-^ = 2h|,(m^^ Mq + M|, -h A\)) — —glM^ — ^ 9i<S 


dt 
MMl 


Svr^-^ = +Mi + Mi + Al) - 3giMi - gfMf - - gfS 




,dMl 

dt 


= 2hl{ml^ +MI + MI + Al) - AglMl + glS 


+ Ml + Ml + Al) - iglMl - g^Ml + i g^S 


( 46 ) 

( 47 ) 

( 48 ) 

( 49 ) 

( 50 ) 

( 51 ) 


Stt^ 


dml[^ 


dt 


3/i|,(m^^ + Mq + Mj) + Ajy) hl{mlj^ Ml + Al) 
-3glMl - glMl glS (52) 


where 

S = ~ + Mq — 2MIj + Ml) — Ml + Ml^ (53) 
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For the bilinear terms in the superpotential we get 


= /i {Shl + 3hl + hl- 3gl - gl) (54) 

= 63 (3hl + K- 3^2 - ^i) (55) 

and for the corresponding soft breaking terms 

dB 

+ Sh'^Ai) + h'^Aj- + 3g‘^M2 + g^Mi (56) 

= 3hljAu+ h‘lAr + 3glM2 +glMi (57) 

The gi are the SU{3) x 57/(2) x 7/(1) gauge couplings and the Mj are the corresponding 
the soft breaking gaugino masses. 


Appendix B: Minimization Procedure 

To minimize the scalar potential we use the procedure developed in refs. P, |^. We 
solve the tadpole equations, eq. (| 5 ), for the soft mass-squared parameters in terms of the 
VEVS and of the other parameters at the weak scale. This is particularly simple because 
those equations are linear in the soft masses squared. To do this we need to know the 
values for the VEVS. These are obtained in the following way: 


1. We start with random values for ht and hf, at Mqut in the range given in eq. (|52D. 
The value of at Mqut is fixed in order to get the correct r mass. 

2. The value of vi is determined from rrib = hbVi/\/2 for nib = 3 GeV (running b mass 
at niz). 

3. The value of V 2 is determined from m* = htV2/\/2 for = 176 ± 5 GeV. If 

vl + vl > v'^ = ^ = (246 GeV)^ (58) 

we go back and choose another starting point. 


4. The value of is then obtained from ^3 



-vf- vl 


We see that the freedom in ht and hb at Mqut can be translated into the freedom in 
the mixing angles (3 and 9. Gomparing, at this point, with the MSSM we have one extra 
parameter 9. We will discuss this in more detail below. In the MSSM we would have 
9 = 7 r/ 2 . 
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After doing this, for each point in parameter space, we solve the extremum equations, 
eq. (^4|), for the soft breaking masses, which we now call {i = Hi, H 2 , L). Then we 
calculate numerically the eigenvalues for the real and imaginary part of the neutral scalar 
mass-squared matrix. If they are all positive, except for the Goldstone boson, the point 
is a good one. If not, we go back to the next random value. After doing this we end up 
with a set of solutions for which: 


1. The Yukawa couplings are determined by the procedure described above. 

2. The other parameters are given by the RGE evolution once the values at Mcut are 
fixed. Notice, however, that these parameters may not satisfy the tadpole equations. 
We will come back to this later. 

3. For a given set of rrif {i = Hi, H 2 , L) each point is also a solution of the minimization 
of the potential. 

4. However, the obtained from the minimization of the potential differ from those 
obtained from the RGE, which we call mj{RGE). 


Our next goal is to find which solutions, for the m? that minimize the effective low- 
energy potential, have the property that they coincide with the m‘f{RGE) obtained, for 
a given unified theory, from the RGE, namely 


ml = ml {RGE) ; i = Hi,H2,L 


(59) 


Following ref. [§] we define a function 


mj 


T] = max 


m: 


'{RGE)' 


ml{RGEy ml 


Vi 


(60) 


Defined in this way it is easy to see that we always have rj >1, the equality being what 
we are looking for. 


We are then all set for a minimization procedure. We want, by varying the parame¬ 
ters at the GUT scale, to get rj as close to 1 as possible. With these conditions we used the 
MINUIT package in order to find the minimum of r]. We considered a point in parameter 
space to be a good solution if < 1 . 001 . 

Before we end this Appendix, let us discuss the counting of free parameters in this 
model and in the minimal N=1 supergravity unihed version of the MSSM. As we explained 
above after requiring the correct masses for the W, t, b and r we get one free parameter 
in the MSSM, tan/3, and two in our model, tan/3 and cos6 or, equivalently, U 3 . As for 
the other parameters we have at the GUT scale one extra parameter, 63 . But we also 
have an extra equation for the tadpoles. So in the end our model has just one more free 
parameter. This can be summarized in the following tables: 
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Parameters 

Conditions 

Free Parameters 

ht, hb, hr, Vi, V2 

mw, rnt, mb, mr 

tan f3 

A, mo, Mi/2, h 

ti = 0,i = l,2 

2 Extra free parameters 

Total = 9 

Total = 6 

Total = 3 


Table 1: Counting of free parameters in minimal N=1 supergravity 


Parameters 

Conditions 

Free Parameters 

ht, hb, hr, Vi, V2, Vo 

mw, rrit, mb, mr 

tan/3, cos6* 

A, mo, Mi/2, F, £3 

t- = 0, i = 1,2,3 

2 Extra free parameters 

Total = 11 

Total = 7 

Total = 4 


Table 2: Counting of free parameters in our model 


Finally, we note that in either case, the sign of the mixing parameter jj, is physical 
and has to be taken into account. 
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